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Abstract 

We prove several inequalities using lowest-order effective field theory for nucleons which give an 
upper bound on the pressure of asymmetric nuclear matter and neutron matter. We prove two 
types of inequalities, one based on convexity and another derived from shifting an auxiliary field. 
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I. INTRODUCTION 



In the effective field theory description of low-energy nuclear matter, nucleons are treated 
as point particles rather than composite objects. While much of the work in the community 
has focused on few-body systems, there has also been recent interest in lattice simulations 
of bulk nuclear matter using effective field theory In parallel with this 

computational effort, effective field theory was also recently used to prove inequalities for 
ihe correlation function of two-nucleon operators in low-energy symmetric nuclear matter 
3]. It was shown that the 5* = 1, / = channel must have the lowest energy and longest 
correlation length in the two-nucleon sector. These results were shown to be valid at nonzero 
density and temperature and could be checked in effective field theory lattice simulations. 
The proof relied on positivity of the Euclidean functional integral measure and is similar 
in spirit to several quantum chromodynamics (QCD) inequalities proved using quark-gluon 
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In this work we prove several new inequalities using effective field theory which give 
an upper bound on the pressure of asymmetric nuclear matter and neutron matter. We 
prove two types of inequalities, one based on convexity and one derived from shifting an 
auxiliary Hubbard-Stratonovich field. We consider two general types of systems, one with 
two fermion species and an SU (2) symmetry and another with four fermion species and an 
SU{2) X SU{2) symmetry. The results we prove are quite general. In addition to nuclear 
and neutron matter, our inequalities apply to systems of cold, dilute gases of fermionic atoms 
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21[ which can be described by the same lowest-order effective field theory. 



II. LOWER BOUND 



Before deriving pressure upper bounds, we first state a general lower bound for the 
pressure. The result is simple and perhaps obvious, but the derivation is useful to help set 
our notation. Consider any system in thermodynamic equilibrium that is invariant under a 
symmetry group S. Let be a symmetric chemical potential which preserves the group S. 
Let /X3 be an asymmetric chemical potential which breaks S and flips sign /is —fi^ under 
some element of S. This means that the pressure P is an even function of fi^. 

Our condition of thermodynamic equilibrium requires that the system is stable and not 
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further separating into regions with more widely different values of /U3. This implies the 
convexity condition, 

> 0. (1) 

Combining this with the fact that P is even in fj,^, we derive the lower bound 

p{|I,lI,)>p{^i,o). (2) 

This lower bound holds for all the systems we consider here. 

III. TWO FERMION STATES - SU{2) 

We consider an effective theory with two species of interacting fermion fields and an 
SU{2) symmetry. Let n be a doublet of fermion fields which we can regard as neutron spin 
states, 



n 



(3) 



i 

We can write the lowest-order Lagrange density in Euclidean space in two equivalent forms, 
jCe = -n[d4 - 2^ + {m% - /i- Atao-s)]^ - \Cnnnn, (4) 

and 

Ce = -n[di - 2^ + {m% - - lJiz(J^)]n - \C'nan ■ nan, (5) 

where 

C = -\C. (6) 

We use a to represent Pauli matrices acting in spin space. // is the symmetric chemical 
potential while nz is the asymmetric chemical potential. We assume the interaction is 
attractive so that 

C < 0, C" > 0. (7) 



A. Two-body operator coefficients 

We can calculate C using a lattice regulator for various lattice spacings, which denote 
as aiattice- For simplicity we take the temporal lattice spacing to be zero. We must sum 
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FIG. 1: Sum over two-particle bubble diagrams. 



all two-particle scattering bubble diagrams, as shown in Fig. ^ and locate the pole in the 
scattering amp litude. We then use Liischer's formula for energy levels in a finite periodic 



box p, |2j, |23| and tune the coefficients to give the physically measured scattering lengths. 
From Liischer's formula there should be a pole in the two-particle scattering amplitude with 
energy 

where ascatt is the scattering length. We can write the sum over bubble diagrams as a 
geometric series. In order to produce a pole at this energy we must have 

1 1 , 1 1 
= lim > f9) 

mjvC ^Tiascatt aiatticeL^ 6 - 2 COS ^ - 2 COS ^ - 2 COS ^ ' 

where aiatuce is the lattice spacing, and the sum is over integer values /ci,A;2,^3 from to 
L — 1. Solving for C gives 

^ / 0.253 \ 

yinascatt aiattice J 

For any chosen temperature and neutron density there is a corresponding maximum value 
for the lattice spacing, aiatuce- The requirements are that the kinetic energy for the highest 
momentum mode must exceed the temperature, and the lattice spacing must be less than 
the interparticle spacing. We therefore have 



hattice ^ lattice) "nr^in — HiaX 



7r- VS^JvT, . (11) 



This sets an upper bound for the absolute value for the scale-dependent coupling C, 

\C\ « I^^Lax : ^ . . • (12) 



I max 

rriN 



47ra. 



This result will be useful for the shifted-field inequalities derived later. 



B. Convexity inequality 



The grand canonical partition function is given by 



Zciti, fis) = / DnDnexp {—Se) = / DnDnexpi / d'^x Ce 



where we use the expression (H)) for Ce, 



= -n[d4 - 2^ + {m% -II- Ii3cr3)]n- - \Cnnnn. 
Using a Hubbard-Stratonovich transformation , we can rewrite Zq as 



Zg oc I DnDnD f exp ( / d^xC^^ 1 , 



where 



^ _ -n[(94 - 2^ + {m% - ^ - lJ.3(J3)]n + Cfnn + \C f . 



Let us define M as the matrix for the part of C^^ bihnear in the neutron field, 



M = - 



^4 - + ^'^^N - - /^3^3) 



Cf. 



We observe that M has the block diagonal form, 



M 



M(/i + /is) 

M(/i - /is) 



where 



+ Cf. 



Since M is real valued, det M must also be real. 
Integrating over the fermion fields gives us 



ZcifJ', fJ-s) cc I DnDnDfexpl I d^xC^ 



D0detM= / DedetM(// + /X3)detM(/i-/X3), 



where DQ is the positive measure 



DQ = Dfexp nC I d^xf 



Using the Cauchy-Schwarz inequality we find 



DQ det M(/i + /ig) det M (/x - ^3) 



< / D0|detM(^ + /i3)detM(/i-/i3)| 



< y y DB [det M(/i + /i3)]y y [det M{fi - /xs)]' 

(22) 



We can now compare the asymmetric partition function to the symmetric partition function 
at chemical potentials yU + ^3 and fx — fj-s, 

ZG{^^, /^s) < V^G(^ + /^3,0)v/^G(^-/i3,0). (23) 

We now use the thermodynamic relation, 

ln^c = S, (24) 

where P is the pressure, V is the volume, and T is the temperature. We find the upper 
bound 

Pifi, fis) < \ [P{l^ + /X3, 0) + P(^ - /i3, 0)] . (25) 
C. Shifted-field inequality 

We start again with the grand canonical partition function 

ZcifJ-, fJ^s) = J DnDnexp {—Se) = J DnDnexp^J (I'^xCe^. (26) 
This time we use the other expression Q for Ce, 



I^E = -n[d4 - 2^ + {m% - /i - fi3cr3)]n - \C'nan ■ nan. (27) 

We can rewrite the grand canonical partition function using three Hubbard- Stratonovich 
fields, 

oc y DnDnD^exp (^j d^x , (28) 

where 



4 = -^[<94 + fi3(T3)]n + iC%- nan - \C% ■ $. (29) 



Let Mo be the neutron matrix without the /iscrs term, 



Mo = - [^4 - 2^ + {m% - /x)J + zC"0 ■ a. (30) 

We note that 

a2Moa2 = Ml, (31) 

where Mq is the complex conjugate of Mg. This means that Mq is either singular, in 
which case det Mq = 0, or has the same eigenvalues as Mq. In all cases det Mq is real. 
Furthermore the fact that a2 is antisymmetric means that the real eigenvalues of Mq are 
doubly degenerate, and so detMp > Q|. 

We now concentrate on the part of £^ that contains ^3 and ^3, 



We can rewrite this as 



where 



— \C' 4>^ + iC (t)'ina^n + ii^na^n. (32) 



-\C'(I)'^ - + iC'<^,na^n + \% (33) 



03 = 03-^^. (34) 

The original contour of integration for ^3 is off the real axis, but we can deform the contour 
onto the real axis. For notational convenience we now drop the prime on ^3 and have 

L% = -n[d4 - + (m% - ^)]n + iC'$- nan - ^C'(j)- $ - ifish + (35) 

The neutron matrix is now Mq, which we have shown has a non-negative determinant. The 
complex phase is contained entirely in the local expression —2/13^3. 
We now have 

Zg oc J DQexp i^J d'^x -«/i303 + 

= ^^P(2^SiT) J DQexp (^-ifi3 j d'^x 03 ) , (36) 
where DQ is the normalized positive measure 



D(P det Mo exp - / d^x V(0) 

DQ = ^ ^ (37) 

/ D(f) det Mo exp ( - / d^x V(0) 
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with 



(38) 



Using pijl we find 



V 



2C" ^ y 



exp( ^^,^^y j J DQexp \^-i^3 J d x (ps 
DQ exp I —ifis / (i'^x 



So we conclude that 



p(^,^3)<^'(m,o) + |^ 



(39) 



(40) 



This upper bound is unusual in that it relates physical observables independent of the cutoff 
scale to the scale-dependent coupling C. By taking the lattice spacing as large as possible, 
we have 

^'=||CLax, (41) 

where |C|j^g^^ was defined in (O, and therefore 



3^1 



2|C|, 



(42) 



As a rough estimate of the quantities involved, we note that for p ~ 0.1 and T < 10 
MeV, \C\^^^ is about 3 fm^. 

As C decreases the upper bound in (pUj) increases. But at the same time the tightness 
of the bound becomes poorer as complex phase oscillations due to the term 



exp 



(43) 



become more significant. The average phase for our functional integral is given by 



(phase) 



DQ exp ( —ips / d^x 03 



exp 



^ lp{p,p,)-P{p,0)-^ 



2C' 



(44) 



Given an estimate of the pressure difference, this relation can be used to predict the 
feasibility of a numerical simulation using this representation of the functional integral. In 
cases where the phase problem is not too severe we can use hybrid Monte Carlo to gener- 
ate Hubbard- Stratonovich field configurations according to the relative probability weight 
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detMo. The phase of the configuration can then be included as an observable using the 
local expression — i/i303. This local expression for the phase could increase algorithmic 
speed by several orders of magnitude. The only known way to compute the phase of matrix 
determinants is LU decomposition, an algorithm which writes a matrix as a product of lower 
and upper triangular matrices. The number of operations for LU decomposition scales as 
A^^, where is the dimension of the matrix. For an lattice the scaling is thus L^^. 



IV. FOUR FERMION STATES - SU{2) x SU{2) 



We now consider an effective theory with four species of interacting fermions and an 
SU (2) X SU (2) symmetry. Let be a quartet of fermion states, which we can regard as 
nucleon fields. 



A^ 



p 




T 


n 







(45) 



We use p{n) to represent protons(neutrons) and t(|) to represent up(down) spins. We use r 
to represent Pauli matrices acting in isospin space and a to represent Pauli matrices acting 
in spin space. We assume exact isospin and spin symmetry in the absence of symmetry- 
breaking chemical potentials, and so the symmetry group is SU{2)j x SU{2)s. 

In the non-relativistic limit and below the threshold for pion production, we can write 
the lowest-order terms in the effective Lagrangian in two equivalent ways. 



-N[d4 



2mjv 



+ {m% - ^)]N - ^Cs{NNf - ICtNBN ■ NaN 



(46) 



or 



-N[di 



2mjv 



+ {m%-^)]N-\C's{NNf 



\C'^NfN ■ NrN 



(47) 



We will introduce symmetry breaking chemical potentials later. We have included both 
three-body and four-body forces. The 5'[/(4)-symmetric three-nucleon force is needed for 
consistent renormalization and has been shown to be the dominant three-body force contri- 



bution 
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With four distinct fermion species there are two irreducible representations of SU{2)j x 
SU{2)s for two fermions in an s-wave, a spin-singlet isospin-triplet (5* = 0) or an isospin- 
singlet spin-triplet (/ = 0). One can show that [7[ 

C[j = —Ct, C'g = Cs — 2Ct- (48) 

In the case of nucleons, one finds that both of the s-wave channels are attractive, with the 
/ = channel being more strongly attractive, 



^scatt '^scatt 



This implies that 

Cs < 3Ct, Ct < 0, (50) 
C's<-C'^, C'^>0. (51) 

For a more general system with four fermion states and an SU(2) x SU{2) symmetry, we 
can interchange the isospin and spin labels so that, without loss of generality, 

' >-4r- (52) 



„/=o - „s=o ■ 

^scatt "'scatt 



In the special case when the scattering lengths are equal, the symmetry group is the full 
Wigner SU{A) symmetry 3^, and the isospin and spin labels can be interchanged. 



A. Two-body operator coefficients 

We determine the two-body operator coefficients in the same manner as before. The only 
difference is that there are now two s-wave channels. The coefficient C in (j^I) is replaced 
by (7-^=° and C^=° where 

C'=' = C's + C'^, (53) 
C^=o = C's- 3C{j. (54) 



We then find 



C's^ ^- + (55) 



C'u ^ 7 V 7 ^ (56) 



O-lattice I \ 47ra.„„tt Hattice. 
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For any chosen temperature and nucleon density there is again a corresponding maximum 
value for the lattice spacing, 



'^lattice ^ ('^/attice )min ~ 



IT 



- V2m^2-, P'^'J • (57) 
This sets a maximum value for the absolute value of the coupling C^, 



\C'u\ « \C'u\ 



-,S=0 



AniN 



lattice J mm 



(58) 



A similar bound for C'g can be made but is not needed in our analysis. 
B. Convexity inequality for /if 

We first consider the case when an asymmetric chemical potential /uf is coupled to the 
nucleon spins. The grand canonical partition function is given by 

Zg = J DNDNexp{-SE) = J DNDNexp (^j d^x , (59) 

where we take the form oi Ce given in (j47|) with an asymmetric spin chemical potential, 

Ce = -N[dA - 2S7 + « - - /^f ^3)]iV - \Cs{NNf - \C'^NfN ■ NtN 

- I^C.iNNf - j,C,iNN)\ (60) 

Using Hubbard-Stratonovich transformations we can rewrite Zq as 

Zgoc j DNDNDfD^exp (^j ct^x C^/^ , (61) 



where 



L{^'^ = -N[d, -^ + {m%-fi- 4as)]N + fNN + iC'J ■ NtN 



2m M 

+ gif)-'^C'J-$. (62) 

In |31[ it was shown that three-body and four-body forces can be introduced without spoiling 
positivity of the functional integral measure. The only requirements are that the three-body 
force is not too strong and the four-body force is not too repulsive. Estimates of the three- 
and four-body forces suggest that these conditions are satisfied. For our analysis here we 
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assume that to be the case, and the function g{f) is a real-valued function which produces 
the two-, three-, and four-body force terms involving NN. 
The nucleon matrix M has the block diagonal structure 



M 



M{fx + nl) 

Mill - fil) 



(63) 



where the upper block is for up spins and the lower block is for down spins. M is a matrix 
in isospin space. 



T. 



We note that 



T2MT2 = M*, 



and so det M > 0. 

Integrating over the fermion fields gives us 



Zaifi^lis) oc J DNDNDfD(j)exp d^x a£ 

= y DO det M = y" De det M(// + //f ) det M(// - //f ) , 



where 



with 



Prom the Cauchy-Schwarz inequality we get 



^g(/u, < \l Zaili + /if > 0)^ Zg{ii - //f > 0). 



We therefore find an upper bound for the pressure. 



(64) 



(65) 



(66) 



(67) 



(68) 



(69) 



P(//, lit) < - [P{iJi + iJil 0) + P{ii - III, 0)] . 



(70) 
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C. Shifted-field inequality for /X3 



We now consider the case with an isospin chemical potential /if. We start with the 
Lagrange density in terms of the Hubbard- Stratonovich fields, 

= -N[d, -£- + {m%-^i- ^iiT,)]N + fNN + zC'J- NrN 

+ g{f) - IC'J ■ (71) 

Let Mo be the nucleon matrix without the //gTs term, 



Mo 

We note that 



+ +f + ^C'J■T. (72) 



T2M0T2 = M;, (73) 



and so det Mq > 0. 

As we did for the two fermion case, we now shift the ^3 field and find the inequality 

F(/i,//f)<F(/i,0) + ^. (74) 
If we take the lattice spacing as large as possible then 

P(^,^,i)<P{^,O) + 0^, (75) 

I ^ I max 

where |C*[;|j^j^^ was defined in (I58|) . As a rough estimate of the quantities involved, we note 
that for p ~ O.Iom and T < 10 MeV, ICrrL ^ is about 0.2 fm^. In this case however the 
situation is complicated by nuclear saturation, and it is not clear that the pionless effective 
theory is applicable. 

V. SUMMARY AND DISCUSSION 

The main results we have shown are as follows. We first considered the two fermion 
system with an attractive interaction and SU{2) symmetry. If p is the symmetric chem- 
ical potential and /is is the asymmetric chemical potential, we proved both the convexity 
inequality 

P(//, 0) < P(//, /is) < \ [P(/i + /is, 0) + P(/i - /is, 0)] , (76) 
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and the shift-field inequahty 

P{fi, 0) < P(/i, fis) < P{fi, 0) + (77) 

We then analyzed the four fermion system with an SU{2)i x SU{2)s symmetry. We 
considered the case when both s-wave channels are attractive and without loss of generality 
assumed the 1 = channel to be more strongly attractive. With /x as the symmetric 
chemical potential and yuf as the asymmetric spin chemical potential we proved the convexity 
inequality 

P(/i, 0) < P(/i, fii) < i [Pi^ + //f , 0) + P(/i - /if, 0)] . (78) 
For non-zero asymmetric isospin chemical potential fiJ^ we proved the shifted-field inequality 

P(//, 0) < P(/i, < Pifi, 0) + (79) 

I ^ I max 

In the Wigner SU{4) symmetry limit, we note that the shift-field inequality ^7^i becomes 
meaningless since |C'^|j^ax ~^ ^- However in this limit we also have the convexity inequality 
for 

P(/i, 0) < P(/i, ^i) < ^ [P{^ + /ii 0) + P(/i - fil 0)] . (80) 

The equation of state for nuclear matter with small isospin asymmetries can be measured 
indirectly in the laboratory by studying nuclear multifragmentation. Of the inequalities 
presented here, the simplest and perhaps most interesting to check is the isospin convexity 
inequality (j8(Jj) in the Wigner symmetry limit. Since much is still unknown about asym- 
metric nuclear matter, this Wigner pressure inequality may be a useful consistency check 
for proposed phenomenological models for asymmetric nuclear matter. 

While some of the inequalities are difficult to observe in nuclear physics experiments, each 
of our results could be tested in the cold Fermi gas system where parameters in the effective 
Lagrangian can be tuned. Such experiments can in principle test the inequalities over a 
range of physical parameters and probe universal results in the limit of infinite scattering 
length and zero range. Although four fermion systems have not yet been produced, these 
may be possible in the near future. 

On the computational side, the inequalities can also be checked by non-perturbative 
lattice simulations. There have been several recent simulations of effective theories on the 
lattice 0,0,0,0]. It will be particularly interesting to look at symmetric and asymmetric 
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nuclear matter in the Wigner symmetry limit, which can be simulated without any sign 
problem. 

It remains to be seen how well many-body nucleon systems can be described without 
explicit pions. Results from [5] for dilute neutron matter suggest that lowest-order effective 
field theory without pions works very well in describing the neutron equation of state. The 
situation for nearly symmetric nuclear matter, however, is less clear due to the effect of 
saturation which requires higher densities. 

With pions included the effective theory action can in general become negative. This 
would in principle invalidate any inequality based on positivity of the action. However it 
has been shown that this sign problem goes away in the static limit js^. Furthermore the 
sign problem has been numerically observed to be small jsj in simulations with neutrons 
and neutral pions for temperatures above 10 MeV and densities at or below normal nuclear 
matter density. If one neglects these sign changes, then the sign-quenched results for the 
effective theory with pions will also satisfy each of the inequalities proven here. 

The author thanks Jiunn-Wei Chen and Thomas Schaefer for several helpful disucssions. 
This work was supported by Department of Energy grant DE-FG02-04ER41335. 
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